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Both in heavy-ion collisions as in magnetars very strong magnetic fields are produced, which has 
f^} its influence on the phases of matter involved. In this paper we investigate the effect of strong 

H ' magnetic fields (B ~ 5mJ/e = 1.7 x 10 19 G) on the chiral symmetry restoring phase transition using 

the Nambu-Jona-Lasinio model. It is observed that the pattern of phase transitions depends on 
the relative magnitude of the magnetic field and the instanton interaction strength. We study two 
specific regimes in the phase diagram, high chemical potential and zero temperature and vice versa, 
which are of relevance for neutron stars and heavy-ion collisions respectively. In order to shed light 
on the behavior of the phase transitions we study the dependence of the minima of the effective 
■ potential on the occupation of Landau levels. We observe a near-degeneracy of multiple minima 

y—t ' with differing occupation numbers, of which some become the global minimum upon changing the 

magnetic field or the chemical potential. These minima differ considerably in the amount of chiral 
i , i symmetry breaking and in some cases also of isospin breaking. 

P-f PACS numbers: 12.39.-x,11.30.Rd,21.65.Qr,25.75.Nq 

4^ I. INTRODUCTION 

£NJ . Recently it has been noted that very strong magnetic fields can be produced in heavy- ion collisions [IHl] • Estimates 
J> ' are that at RHIC magnetic fields are created of magnitude 5.3m^/e = 1.8 x 10 19 G and at LHC of 6m^/e — 2x 10 19 G, 
and there are even higher estimates [||. Also, certain neutron stars called magnetars exhibit strong magnetic fields, 
between 10 14 — 10 15 G 0, @. These fields occur at the surface, probably in the much denser interior even higher fields 
are present. Using the virial theorem it can be derived that the maximal strength is 10 18 — 10 19 G 7]. If one assumes 
• ■ that the star is bound by the strong interaction instead of by gravitation, this limit can be even higher. 

In both neutron stars and in heavy-ion collisions it is expected that quark matter plays a role. Therefore it is 
interesting to study how this form of matter behaves in a strong magnetic field. Two different regions in the QCD 
phase diagram are of relevance here. Heavy-ion collisions probe the low chemical potential and high temperature 
regime, for neutron stars it is the other way around. In this paper the effect of very strong magnetic fields will be 
investigated in both regimes. 

Much work has been done on how an external magnetic field changes nuclear matter, for a review see Ref. [|. The 
behavior of ordinary quark matter has been studied using the Nambu-Jona-Lasinio (NJL) model, see for example 
Refs. [9l-[l^| and recently also in the linear sigma model coupled to quarks [lj|. Most studies investigate the one 
and two flavor cases, but recently also the three- flavor case has been investigated [2(J[2l|]. At high quark chemical 
potential, it is believed that the ground state is a color superconducting phase. The effects of an external magnetic 
field on such a phase are discussed in Refs. [23 - l28j . Here color superconductivity will not be considered. 

In this paper we study the chiral symmetry restoring phase transition, which is strongly influenced by an external 
magnetic field. From studies in the NJL model it is known that a magnetic field enhances the chiral symmetry 
breaking [9J, this is related to the phenomenon of magnetic catalysis of chiral symmetry breaking, introduced in 
Refs. [29 31], further studied for the NJL model in Refs. [lMi and for QED in e.g. Refs. [32l - [36| . where also the 
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generation of an anomalous magnetic moment was pointed out [35L [36} . This enhancement of chiral symmetry breaking 
can be understood as follows, the B-ficld anti-aligns the helicities of the quarks and antiquarks, which are then more 
strongly bound by the N JL-interaction Q . The phenomenon of magnetic catalysis of chiral symmetry breaking leads 
to interesting behavior, since it allows for phases with broken chiral symmetry and nonzero nuclear density for a range 
of chemical potentials and magnetic fields [lil [lH Kl\ . In such a phase nonperiodic magnetic oscillations occur, which 
means that the constituent quark masses are strongly dependent on the magnetic field, and consequently also other 
thermodynamic parameters. 

In all studies of the influence of magnetic fields on chiral symmetry breaking up to now, the effects of instantons 
have not been studied explicitly, i.e. as a function of instanton interaction strength. Magnetic fields and instantons 
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can lead to combined effects. In Ref. [3j it is shown that variations in topological charge, which induce variations of 
net chirality, in a strong magnetic field gives rise to an electrical current. This effect is known as the chiral magnetic 
effect and could perhaps be observed in heavy-ion collisions. Variations of topological charge can for instance be 
created by instantons. 

Here a related study will be performed. We will investigate the combined effect of instantons and a strong magnetic 
field on quark matter using the NJL model. In this model instantons induce an extra interaction, the 't Hooft 
determinant interaction, which leads to a mixing between the different quark flavors. Following the analysis of 
Refs. [371 . |38[, the strength of the instanton interaction is set by the dimensionless parameter c. For c = there is 
no contribution and the quarks are fully independent. Because of the difference in charge of the quarks, the phase 
transitions are decoupled. The other extreme case is c = 1/2, which is actually the most studied case. The quarks 
are then fully mixed, the constituent quark masses are equal and the phase transitions will always coincide. In Ref. 
[l4j this case is studied in the chiral limit. It is observed that for a range of typical value for the coupling constant, 
phases with broken chiral symmetry and nonzero nuclear density arise. 

In general there is a competition between the magnetic field, which tends to differentiate the constituent quark 
masses for different flavors, and the instanton interaction which favors equal constituent quark masses. In this work 
this competition is studied. Apart from studying the ground state as function of magnetic field and chemical potential 
for various characteristic values of c, we also look at the local minima of the effective potential and the corresponding 
occupation of Landau levels. It is found that in the neighborhood of the chiral phase transition the phase diagram 
develops metastable phases, differing in the number of filled Landau levels. Some of these local minima become the 
global one upon increasing the magnetic field or chemical potential, but not all of them do. These phases can have 
rather different values for the constituent quark masses, in other words, display significantly different amounts of 
chiral symmetry breaking. Unlike in the case of c = 1/2 which is isospin symmetric, in these phases the values of the 
two constituent quark masses can be very distinct, which corresponds to large isospin violation. Furthermore, we find 
that for a realistic choice of parameters, appearance of phases of broken chiral symmetry and nonzero nuclear density 
requires not too large instanton interaction strength, i.e. c<0.1. 

As mentioned, we also investigate the role of nonzero temperature at zero chemical potential, which is of relevance 
for heavy-ion collisions. Without magnetic field the chiral symmetry restoring phase transition at finite temperature 
is a crossover. In the linear sigma model coupled to quarks it has been observed that the magnetic field turns it into 
a first order transition [l!| . We will see that this is not the case for the NJL model. 

This paper is organized as follows. First we derive the effective potential of the NJL model in the mean-field 
approximation in a magnetic background. Then we discuss the phase diagram as a function of chemical potential, 
concentrating on the phase with nonzero nuclear density and chiral symmetry breaking. We continue with discussing 
the temperature dependence and end with conclusions. 



II. EFFECTIVE POTENTIAL OF THE NJL MODEL WITH A MAGNETIC FIELD 

In this section the effective potential of the NJL model in a strong magnetic background is derived, this analysis 
is based on Refs. [H, [l9| . The effective potential derived in this way is equal to the one in Refs. [TUT! using the 
Fock-Schwinger proper time method. 

First to set our notation, we briefly review the NJL model introduced in Refs. [39|, |4(J. It is a quark model with 
four-point interactions, the gluons are 'integrated out'. In this paper the following form of the NJL model is used, in 
the notation of Ref. [38| 

£njl = $ (*7 m ^m -m- /ryo) ip + jC Sq + C dc t, (1) 

where m = m u = nid is the current quark mass and /i = \x u = fid is the quark chemical potential. Note that in 
contrast to Ref. (38j . here the current quark masses and both chemical potentials are taken equal. Furthermore, 

c qq = d [(^KVO 2 + OKz 75l /o 2 ] , (2) 

is the attractive part of the qq channel of the Fierz transformed color current-current interaction [4l| and 

Act = 8G 2 det (M L ) + h.c, (3) 

is the 't Hooft determinant interaction which describes the effects of instantons [12, EH . Note that this interaction is 
flavor mixing. In the literature G\ and Gi are often taken equal, which means that the low energy spectrum consists 
of a and 7r fields only, but here we will allow them to be different. We will restrict to the two flavor case, using r a 
with a = 0, 3 as generators of U(2). 
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The symmetry structure of the NJL model is very similar to that of QCD. In the absence of quark masses 
and the instanton interaction, there is a global SU(3) c xU(2)iXU(2)^j-symmetry The instanton interaction 
breaks it to SU(3) c xSU(2)j / xSU(2)ijxU(l)B- For nonzero, but equal quark masses this symmetry is reduced to 
SU(3) c xSU(2)y xU(1)b. If a magnetic field is turned on, the symmetry is reduced to SU(3) c xU(l) 2 due to the 
differences in charge. 

We choose the parameters as in Refs. [13, This means we write 



G 1 =(l-c)G , G 2 =cG , (4) 

where the parameter c controls the instanton interaction strength, while the value for the quark condensate (which 
is determined by the combination G\ + G2) is kept fixed. For our numerical studies we will use the following values 
for the parameters: to = 6 MeV, a three-dimensional momentum UV cut-off A = 590 MeV/c and GqK 2 = 2.435. 
These values lead to a pion mass of 140.2 MeV, a pion decay constant of 92.6 MeV and finally, a quark condensate 
(uu) — (dd) = (—241.5 MeV) 3 [l?]], all in reasonable agreement with experimental determinations. 

To calculate the ground state of the theory, the effective potential has to be minimized. In this section the effective 
potential is calculated in the mean-field approximation. We will assume that only the charge neutral condensates 
(ipTQtp) and (ipT3ip) can become nonzero. 

To obtain the effective potential in the mean-field approximation, first the interaction terms are "linearized" in 
the presence of the (ipTQip) and (V'73'0} condensates (this is equivalent to the procedure with a Hubbard-Stratonovich 
transformation used in Ref. 38]) 



leading to 



(lpT a lp) 2 ~ 2 (ipTalp) 1pT a lp - (lpT a lp) 2 , (5) 



£ NJL = * (n - M - M7o) * - ^gf^ - w ^ &; , (6) 



where M = M tq + M 3 t 3 and 

M = to - 2G (iM,V) , 

M 3 = -2(1 - 2c)G (iprsip) . (7) 

Now the Lagrangian is quadratic in the quark fields, so we can integrate over them. After going to imaginary time, 
this results in the following effective potential in the mean-field approximation 44] 

V ' f=u P0 = (2n+l)TrT J V ; 

where we have introduced the constituent quark masses M u = Mq + M 3 and Md = Mq — M 3 . 

As we have mentioned earlier, often Gi is taken equal to G2, which is the case c = 1/2. This choice implies that 
M3 is then always equal to 0, i.e., the constituent quark masses are equal. Note that the reverse is not true. If the 
assumption is made that the constituent quark masses are equal, c = 1/2 or {t]-!T 3 ip) = or both. However, for M 3 = 
changes in c cannot be noticed because only the combination Gi + G2 occurs. Hence, the conclusion is that (^T 3 tp) = 
and c can be any value. Instanton effect are nevertheless present (G2 can be nonzero after all) and 77 and ao mesons 
can still be present in the spectrum. The ratio G2/G1 can simply not be determined if A/3 = 0. Since magnetic fields 
affect the two flavors differently because of the difference in charge, isospin breaking effects are expected and it is 
unnatural to choose M 3 = from the start or to take c = 1/2. Using the strange quark condensate it was argued in 
Ref. [3?! that a realistic value of c would be around 0.2. 

A. Including a magnetic field 

Now we include a magnetic field, which changes the dispersion relation for the quarks in the following way 

p 2 Qn =p 2 z +M 2 + (2n+l-a)\q\B, (9) 

where n is the quantum number labelling the discrete orbits, a the spin of the quark, and q is its charge. The integral 
over the three-momentum is modified as 

d 3 p ^ \q\B ^ f dp 



(2tt) 3 2tt ^ J 2tt 
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Performing the sum over the Matsubara frequencies gives the following effective potential [14j, |l_5j, |l£ 



V 



(M - mf 



Ml 



4G 



4(1 - 2c)G 2tt 



-Tin 



l + e -[E pJ (B)- f A/T | ; 



(11) 



where E p j(B) = \jV 2 + (2n + 1 — <j)\qj\B + Mj. Following the analysis of Ref. [18J this potential can be split in 

three pieces, a part that is independent of external parameters, a part that only depends on the magnetic field and a 
part that depends on the magnetic field, chemical potential and temperature 



with 



Vn = 



V = V + V 1 (B)+V 2 (B, M ,T), 
(M -m) 2 , 



(12) 



Ml 

AG T 4(1 - 2c)G 

d 



2tt 2 



c'(-i,*/)-2<4- a, /) lni / + x 



1 + e -[^ P ,/(S)+H/T 



+Tln [l + e -[^,/(B)-M]/rl j ; 



(13) 
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where we have defined Xf — 2\q/\B an< ^ C(~^-> x f) — ^dz^ \z=—i with C,(z,xj) the Hurwitz zeta function. We have 
neglected xj independent terms in Vi(-B) (including a UV divergent one). 

The term Vo is divergent and needs to be regularized. Here a conventional three-momentum UV cut-off is used, 
yielding the expression 



Vo 



(Mo - mf 
4G 



Ml 



4(i-2c)g 8ti - 2 j^j tl y f yM f 



A 2 
M 2 



A {Mj + 2 A 2 ) 



I A 2 
M 2 



(14) 



The expression 1, a;/) in Vi(B) can be written in a more convenient form by differentiating and integrating the 
function with respect to Xf\ 



C'(-l, x f ) = C'(-l, 0) + ^ - ^ - ^ log(2vr) + ^- 2 > (*/), 



(15) 



where ip( m \xf) is the m-th polygamma function. The term £'(—1,0) is independent of xj and will therefore not be 
taken into account. The remaining expression is amenable to numerical evaluation. 
The summation over a and n in V2(B, n, T) can be rewritten as 



V 2 (B, / x,T) = -^^(2-4o)k/|B / ^{nn 



l + e 



-[£ p , fc (T)+ M ]/T 



Tin 



l + e -[^(T)-M]/rl \ (16) 



where = yp 2 + M 2 + 2k\qf\B and fc denotes the Landau level, which has degeneracy (2 — 5ko). 
At zero temperature, V2 can be simplified to 

N r 
'27 



V 2 (B, M ,0) = -^£(2-4o) / ^(Ai-^k)^-^] 



£ J] (2-M^(M-s/(fc,S)) 



/=« fc=0 



Igfj-ggg 
4tt 2 



nJ^ 2 -s 2 (k,B)-s 2 f (k,B)\n 



H+^Jfi 2 -sj{k,B) 
s/(fc,-B) 



(17) 
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where Sf(k, B) = \ M"i + 2\qf\Bk and kf^ max is the upper Landau level, defined as 



2\q f \B 



(18) 



where the brackets indicate the floor of the enclosed quantity. 

We will now use these expressions in a numerical study of the minima of the effective potential, performed along 
the lines discussed in Refs. [38, Hi. 



III. RESULTS 



We start with considering the case of (i = 0, T = 0, and c = 0. Fig. Q] shows the results for this unmixed case. 
The magnetic field enhances M u and M^, which are proportional to {uu} and {dd}, respectively, consequently the 
chiral symmetry breaking is enhanced |9|. Because of the charge difference of the quarks, the -B-dependence of the 
constituent quark masses is not equal. Nonzero c will cause mixing and will bring the masses closer to each other. As 
discussed, at c = 1/2 the constituent quark masses are exactly equal. 




FIG. 1: The dependence of the constituent quark masses M u and M d on the magnetic field B. 



A. Nonzero chemical potential 



In this section we turn to the phase structure near the phase transition at nonzero chemical potential and zero 
temperature. From Refs. pril - l49l | it is known that when the isospin chemical potential is nonzero, it is possible to have 
two phase transitions at low temperature and high baryon chemical potential. Here we study a similar case, instead 
of nonzero isospin chemical potential, we allow for nonzero magnetic field, we will see that also here the possibility of 
separate phase transitions for the two quarks arises. We will take equal chemical potentials for the quarks, but the 
magnetic field acts effectively like a nonzero isospin chemical potential due to the difference in charge of the quarks. 
Instantons cause mixing between the quarks, if the mixing is strong enough, the two separate phase transitions merge 
into one. This was extensively investigated in Ref. [46j for the nonzero isospin chemical potential case. 

From Ref. 0[ll[l7|], where the NJL model in the chiral limit was studied, it is known that Landau quantization 
induces a more complex phase structure. Apart from the usual phase of broken chiral symmetry with zero nuclear 
density, there is also the possibility of such a phase with nonzero nuclear density. Here we perform a more detailed 
study of this case, which is a characteristic phenomenon at nonzero chemical potential and sufficiently strong magnetic 
fields (cf. Eq. (UU) below). 
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1. The c = case 

When the determinant interaction is turned off, the up and down quarks are decoupled. This leads to the possibility 
of separate phase transitions for the quarks. In Figs. [5] and |3] we show the constituent quark mass of the up and down 
quark respectively as a function of quark chemical potential and magnetic field. As expected, the two quarks have 
decoupled behavior. 

Let us first discuss the behavior of the up quark. At low chemical potential we have the "standard" chiral symmetry 
breaking NJL ground state with empty Landau levels (LL). Following the nomenclature of Refs. [H], EE] where the 
c = 1/2 case was studied in the chiral limit, this is called phase B. Note that this phase always has zero nuclear 
density. At high chemical potential chiral symmetry is restored, up to the explicit breaking. In this approximate 
symmetric phase magnetic oscillations can be seen in the constituent quark masses, caused by Landau quantization. 
These oscillating phases are denoted by Ai, where i gives the number of filled LL. As these phases have occupied LL, 
they have nonzero nuclear density. The nuclear density of level k is given by [l8j 

p ftk (B, m) = (2 - SmWO* s f (k, B)) l -V^^i-s}(k,B). (19) 

In the chiral limit the constituent quark masses vanish in the Ai phases. 

The oscillations are due to the de Haas- van Alphen effect, which in QED and in the two-flavor NJL model for 
c = 1/2 in the chiral limit lead to second order transitions between the Ai phases However, with our choice of 
parameters the transitions are weakly first order. In the chiral limit they become second order, like for c = 1/2, as 
can for instance be seen in the nuclear density. For completeness, we mention that in the color superconducting case 
of Ref. [27, 28] the oscillations in the gap parameter are seen to be continuous, but second order transitions can occur 
when neutrality conditions are imposed. 

For B larger than 4.5m^/e an interesting intermediate phase arises, where the up-quark jumps as a function of \x 
first to a phase with a still rather large constituent mass and then to phase A\. This intermediate phase is called Co 
in the language of Refs. [3, EH] and corresponds to a phase of broken chiral symmetry having nonzero nuclear density 
and a filled zeroth LL. So the essential difference between this phase and Aq is the breaking of chiral symmetry. For 
smaller values of the coupling constant Go also the phases Cfc with k > (which are similar to Co but with more 
occupied LL) occur. The transitions between the Cfc are first order, furthermore, they are nonperiodic in the sense 
that the difference between the transitions is -B-dependent as the constituent mass strongly depends on B [3]. If 
we are in this phase Co and increase the magnetic field, the constituent quark mass decreases, eventually becoming 
almost zero, this can be interpreted as a crossover to Aq. In the chiral limit the crossover becomes a second order 
transition. Finally, we would like to note that already at B = 4mJ/e the phase Co exists as a metastable phase (we 
will discuss this in more detail later). 

The qualitative behavior of the down quark is very similar, as the quarks only differ in charge, consequently Fig. [3] 
can be directly obtained from Fig. [5] by multiplying B by 2, for ease of comparison we show both figures. If one 
compares the two figures, one can immediately see that there are large regions where the constituent quark masses 
are considerably different. This is equivalent to a large nonzero (tpTsip) condensate, i.e., spontaneous isospin breaking. 
This will influence the behavior of the mesons accordingly, for example the masses. 

Eventually, if one keeps increasing the magnetic field, the phase transitions of the quarks will take place at (almost) 
the same chemical potential and there will be no spontaneous isospin breaking. 

2. Thec^O case 

In this section the consequences of the instanton interaction is studied, i.e., the parameter c is varied. Increasing c 
will cause mixing between the constituent quarks, which tends to bring the constituent quark masses together. Around 
the phase transition there is a competition between the effect of the magnetic field and the instanton interaction. 

The competition is illustrated in Fig. [4j where the constituent quark masses are plotted as a function of the quark 
chemical potential for three characteristic values for c, c = 0, 0.03, and 0.1 with B = 5m^/e. The qualitative behavior 
for different values of the magnetic field is similar. One can see that when c ^ 0, the phase transitions are indeed 
coupled. Furthermore, one observes that the two phase transitions merge into one when c is increased and that the 
phase Cp disappears. Qualitatively the behavior is similar to the case of nonzero isospin chemical potential studied 
in Ref. (46|, but in that case the phase Co does not exist. 

When the coupling constant Go is lowered, it is possible to have Cfc phases at c = 1/2, as in Ref. [l4j]- Compared 
to the chiral limit studied there, the region of the phase diagram with C* phases increases for m ^ 0. 

More insight into the phase structure and phase transitions is obtained by looking at the behavior of local minima 
of the effective potential. Near the phase transition at these (large) magnetic fields, metastable phases arise. These 
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FIG. 2: The dependence of the constituent up-quark mass on B and /j,. Ai, B and Co denote the different phases using the 
scheme of Refs. 0, [3]. Chiral symmetry is broken in phases B and Co, phases Aq and Co have nonzero nuclear density. 




FIG. 3: Same as Fig. [2] now for the down quark. 

phases differ in the number of filled LL. Let us take as an example the c— 1/2 case, which is the easiest to discuss, 
as the effective potential is then only a function of M u = — M . In Fig. [5] we show the effective potential as a 
function of M with /i = 378 MeV and B = 5m^/e. At these values four minima can be seen, the global minimum is 
the phase in which the chiral symmetry breaking is largest, i.e. minimum 4. When /i is increased, minimum 1 will 
take over, which is A\ for the up quark and Ai for the down quark. The other two local minima never become the 
global one for our choice of Go, but as they are almost degenerate with the other minima (also for other values of c), 
they are nevertheless important. These local minima correspond to CV- phases and can become the global minimum 
when Go is lowered. 

Similar results hold for c ^ 1/2, also then metastable phases exists with different fillings of LL. In this cases some 
of the Ck phases can become the global minimum, as we have seen for c = 0. Like before, the number of such states 
depends on the choice of the other parameters. 
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FIG. 4: The dependence of the constituent quark masses on the quark chemical potential for B 
Solid lines denote the up quarks, the dashed lines the down quark. 



5ra^/e and various c values. 



As the metastable phases differ in the values of the (uu) and (dd) condensates at small c, they again represent rather 
large broken isospin and will lead to different meson masses. Whenever the system is passing the phase transition, it 
could be trapped in one of those metastable phases for some time and consequences from the changing meson masses 
can arise, for example, enhancing or suppressing certain decays. 



B. Nonzero temperature 



In this section the temperature dependence of the ground state is investigated at zero chemical potential, but with 
a magnetic field. As the instanton interaction does not influence the temperature dependence much, we only consider 
c = 1/2 for simplicity. Ref. [l9| found in the linear sigma model coupled to quarks, the usual crossover becomes a 
first order transition at very high magnetic fields. However, we find that this is not the case in the NJL model. 

In Ref. [l^| only the lowest Landau level was taken into account. Here more Landau levels are included, so the 
effect of the higher Landau levels can be investigated in the NJL model. Since the levels with large k are exponentially 
suppressed, the summation can be truncated in Eq. (I16[) , we will denote the largest k with fc trunc . The value of k tmnc 
depends on the temperature, constituent quark mass, chemical potential and magnetic field considered. If M and T 
are increased or if B is decreased, fc tr unc has to be increased. 

In Fig. IS] we show the temperature dependence of the constituent quark mass at B = 15m^./e for four different 
values of fctrunc- The 13 levels case is chosen such that the error is less than 1 percent at M = 450 MeV, T = 450 MeV. 
From the figure it can be inferred that taking more Landau levels into account, makes the crossover sharper. Also, 
there is a significant difference between including the zeroth and first Landau level. It is clear that including more 
Landau levels, influences the details of the transition. However, the qualitative aspects of the phase transition are not 

-0.999 
-0.9992 
-0.9994 
-0.9996 
_ -0.9998 

> -1 

> -1.0002 
-1.0004 
-1.0006 
-1.0008 

-1.001 

50 100 150 200 250 300 350 400 450 
M (MeV) 

FIG. 5: The normalized effective potential at the values B = 5mJ/e, fj, = 378 MeV and c = 1/2. There are four minima. The 
numbers below the minima denote the LL occupied for each quark. Note that the minima are almost degenerate. 




u: none 
d: none 
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changed. 
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FIG. 6: The temperature dependence of the constituent quark mass for strong magnetic held (B = 15m^/e) and various fc trU nc 
values. 

In Fig. [7] the temperature dependence of the constituent quark mass for different values of the external magnetic 
held is shown. The phase transition remains a crossover, in contrast to the results in the linear sigma model coupled 
to quarks. This difference is important, as a hrst order phase transition allows for meta-stable states, whereas a 
crossover does not. 

In Fig. [3] the results in the chiral limit are shown, where the transition remains a second order phase transition, 
like it is the case at zero magnetic held and conhrms the results of [TtJ who calculated the phase diagram in a strong 
magnetic held in the chiral limit using the Fock-Schwinger proper time method. Note that the critical temperature 
increases slightly with increasing magnetic held. 




50 100 150 200 250 300 350 

T (MeV) 



FIG. 7: The temperature dependence of the constituent quark mass for various strong magnetic fields. We also indicate fctn 



IV. CONCLUSIONS 

The effect of a strong magnetic held on quark matter has been investigated in the NJL model in two regimes, 
zero temperature and finite chemical potential and vice versa. The hrst regime is of relevance for (the interior of) 
magnetars and the second for heavy-ion collisions. 

At very high magnetic fields, when M ~ 2\q\B ss fi, the phase structure shows a variety of phases and phase transi- 
tions due to Landau quantization. As a function of chemical potential, more phase transitions occur, corresponding to 
Landau levels filling up successively. Due to the difference in charge, this pattern is different for the two quark flavors. 
When there is no mixing in the absence of the instanton interaction, the two patterns are uncoupled. This generally 
leads to rather different constituent quark masses, or equivalently, spontaneous isospin breaking (ipr^ip) ^ 0. This 
affects the mesons inside the medium, for example their masses. It was found that for a realistic choice of parameters 
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FIG. 8: Same as Fig. now in the chiral limit. 

in the NJL model such a phase of broken chiral and isospin symmetry arises around B — 4.5m^/e, but it is already 
present as a metastable phase for lower magnetic fields. 

When the instanton interaction is included, a competition occurs between the strength c of this interaction and the 
magnetic field. This reduces the region in the phase diagram with large {if>T3ip}. For c sufficiently large it disappears 
entirely, leaving only one phase transition. However, around this transition the phase structure is still rather complex 
regarding metastable phases, which are characterized by different fillings of Landau levels and which differ only slightly 
in energy, but much in the amount of chiral symmetry breaking. For lower values of c some of these near-degenerate 
minima can also differ considerably in the amount of isospin breaking. 

Finally the role of temperature was studied at zero chemical potential. In Ref. [19j it was found in the linear sigma 
model coupled to quarks that a strong magnetic field changes the usual crossover as a function of temperature into 
a first order transition. In the NJL model it was found that the crossover remains a crossover. Also it was found 
that including higher Landau levels in the calculation of the effective potential changes the details of the crossover, 
it becomes sharper, albeit the qualitative aspects of the transition are not changed. The difference between the two 
models is important, as the first order transition allows for metastable phases, while a crossover does not. 
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